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1 Abstract 

Flame structure and ignition and extinction processes associated with a strained fuel 
strip are explored numerically using detailed transport and complex kinetics for a propane- 
air reaction. Ignition modes are identified that are similar to those predicted by one-step 
activation energy asymptotics, i.e., modes in which diffusion flames can ignite as independent 
or dependent interfaces, and modes in which single premixed or partially premixed flames 
ignite. These ignition modes are found to be dependent on critical combinations of strain 
rate, fuel strip thickness, and initial reactant temperatures. Extinction in this configuration 
is seen to occur due to fuel consumption by adjacent flames, although viscosity is seen to 
have the effect of delaying extinction by reducing the effective strain rate and velocity field 
experienced by the flames. 


2 Introduction 

Ignition and extinction processes associated with a single strained fuel-air interface have 
been examined extensively over the years using single step activation energy asvmptotics[l- 
4] as well as reduced kinetics simulation[5-7] and full combustion chemistry simulation [8, 
9]. These processes may be understood in terms of the classical “S-shaped curve” [2, 3, 10], 
representing the thermal runaway leading to ignition as strain rate is relaxed and flame 
extinction as strain rate is increased. The concept of ignition delay due to straining has 
been applied to a variety of situations, including an understanding of flame liftoff[4, 11] as 
well as ignition delay processes in jets, shear layers, and vortical structures[6, 12-14], 

In many complex flow situations, strained fuel-air interfaces form[15] and are brought 
closely together into layer-like turbulent flow structures which may be viewed as alternat- 
ing fuel-oxidizer strips undergoing local strain. These dissipation and reaction layers are 
comprised of laminar “flamelets” which are used to model turbulent non-premixed combus- 
tion[16]; it has been suggested that the existence of these layers is independent of the extent 
of chemical non-equilibrium in the flow[17]. 

While extinction of the interfaces bounding finite regions of fuel has been considered 
in the fast chemistry limit[18, 19] as well as using finite rate single step kinetics[20] due 
to fuel burnout, the interesting problem of ignition delay associated with a strained fuel 
layer, bounded by oxidizing species, has only recently been studied via single step activation 
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energy asymptotics[21]. This study examined how the ignition process is affected as strained 
interfaces are brought closely together due to straining. Three possible modes of ignition are 
determined: one in which the two interfaces ignite independently to form diffusion flames, 
one in which the two interfaces ignite in a dependent fashion and in which ignition ma_v be 
actually augmented by strain due to enhanced thermal feedback, and one in which ignition 
occurs to form a single, premixed flame at very high strain rates before ignition is completely 
prevented. In contrast to a single, isolated interface in which ignition can be prevented by 
over-matching heat production with heat convection due to strain, ignition of a strained fuel 
strip can also be prevented if the finite extent of fuel is diluted by oxidizer more quickly 
than heat production can cause a positive feedback thermal runaway. These behaviors are 
dependent on the relative sizes of time scales associated with species and heat diffusion, 
with convection due to strain, and with the chemical reaction. Critical strain rates leading 
to complete ignition delav can be considerably smaller for the fuel strip than those for single 
interfaces as the fuel strip is made thin in comparison to diffusion and chemical length scales. 

When critical strain rates leading to ignition prevention for two adjacent species interfaces 
are smaller than those required to delay ignition of a single isolated interface, certain benefits 
may be derived. If ignition is delayed due to strain in a turbulent reactive flow, fuel and air 
can experience substantial mixing, so that when the strain field relaxes, ignition maj occur in 
a premixed mode, potentially reducing NOx production that would occur in stoichiometric 
diffusion flames[22]). The strained fuel strip is in fact a component of a low NOx, lobed 
fuel injector currently under examination[23, 24], Streamwise vorticitv generated by the 
injector creates a flowfield of high strain, acting to delay the ignition of injected strips of 
fuel into a coflowing airstream. These issues are what motivated the earlier study of ignition 
delay in strained fuel layers via activation energy asymptotics[21]; they continue to motivate 
the present, more quantitative study which examines ignition, burning, and extinction in a 
propane fuel strip surrounded by air using a full kinetic mechanism. 


3 Problem Formulation and Solution Procedure 

The geometry of the strained fuel strip problem is shown in Figure 1. The strain rate 
imposed in the problem is formulated through boundary conditions as y — ► +oo, that is, 
u(x,y -> +oc ,t) = ex and v{x,y -> +oo,f) = -ey. It should be noted that the boundary 
conditions as y — > — oo are the same as those a y > +oo due to the presence of air in 
both limits, in contrast to the classical counterflow diffusion flame problem with differing 
initial densities of fuel and air. The strained fuel strip problem is actually richer from a 
phenomenological point of view than is ignition of the single diffusion flame in that additional 

modes of ignition are possible, as noted in [21]. 

The governing equations for the strained fuel strip problem (mass, momentum, energy, 
and species conservation) may be solved in general in the present scheme in conservative 
form: 
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( 1 ) 


Ut + [F(U)] x + [G(U)] y = [£,(#)] r + [G v (t/)] y + S 

where the vector U (containing the density p, velocity components u and u, energy per 
unit volume e, pressure p, and mass fractions Yi), the convective flux vectors F and G, the 
diffusive flux vectors F v and G v , and the source vector 5 are defined as follows: 
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Mass production rate terms w* are evaluated in (4) in standard fashion, assuming an Arrhenius- 
tvpe reaction rate coefficient dependence on temperature. Inviscid as well as viscous forms 
of the equations have been solved in one and two dimensions, respectively[25]. 

If a one-dimensional formulation is made and if pressure is assumed spatially invariant, 
as done by Darabiha[26] and by Petrov and Ghoniem[27], for example, with the simplified 
boundary conditions, the problem essentially reduces to a set of energy and species conser- 
vation equations: 


u t + [G{uj\ = + [£,(#)] y + s + c 


( 5 ) 
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where the heat flux component, neglecting radiative heat transfer, is given by: 


( 6 ) 


(7) 


NS 

Qy = \T y + pY,hkD km (Y k ) y ( 8 ) 

t=i 

The above equations may be solved numerically. The energy conservation equation is 
solved in terms of the energy per unit volume e (or actually, the enthalpy per unit mass, h). 
This is in contrast to the common assumption of a “locally” calorically perfect gas where 
h = c p T is assumed in the governing equations, employing table look-ups to obtain the 
dependence of c p ,; for each species on temperature in the governing equation and averaging 
to obtain c p . While for many combustion problems the usual approach yields accurate results, 
it has been recently shown [29] that such assumptions in the vicinity of sharp gradients such 
as flames can yield errors in wave speeds and hence in the solution of the problem. The 
present approach avoids this difficulty. 

Formulation of the ratio of specific heats for the mixture of gases present in a given 
computational cell is done so as to avoid thermodynamic inconsistencies at contact discon- 
tinuities. These inconsistencies are known to generate temperature spikes and numerical 
oscillations[30-32]. In the present procedure we use the relation 


7 = 1 + 



(9) 


for the ratio of specific heats of the gas mixture in a given cell; this assumes that all species 
Y k with molecular weights W k are at the same temperature within the cell. Although this 
assumption is not as accurate as evolving 7 from partial pressures of the species, as done by 
Ton [30], experience has shown that if a sharp contact surface rapidly fades away due to the 
presence of highly diffusive species, then using equation (9) will result in confined spikes in 
temperature for only very short periods of time. 
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The method of Strang splitting or operator splitting via fractional steps[33] is employed 
to solve the system of equations (1) or (5). This procedure enables the system to be split 
into two separate equations, one which only includes the advection-diffusion terms and one 
which only includes the reaction rate source terms, respectively: 


U t + 


F(U)} X + [G(U)] = [Fv{U)] x + [Gv(U)] y 


( 10 ) 


U t = S ( 11 ) 

An equivalent form of equations (10) and (11) may be written for the inviscid formuation 
of the problem, equation (5). Each relation above represents a set of equations that has 
been studied extensively. The present method uses a first-order upwind scheme[34] to solve 
the inviscid form of the advection-diffusion equation (10) and the second-order McCormack 
predictor-corrector scheme[36] to solve (10) in the full viscous form. It employs a stiff ODE 
solver (VODE[35], an improved version of LSODE(37]) accompanied by the CHEMKIN-II 
package [38] to solve equation (11). 

The above formulation has been demonstrated to be relatively efficient as well as accurate 
in its prediction of the structure of strained diffusion flames and fuel strips. The present 
calculation employs the propane-air chemical mechanism of [39] to study ignition, burning, 
and extinction of a strained propane fuel strip. 


4 Results 

Results for the ignition, burning, and extinction processes associated with the propane 
strip burning in air are shown in Figures 2-6 for the inviscid case. Figure 2 is a plot of 
temperature as a function of distance across the fuel strip and as a function of time, for a 
fuel strip (T = 400° K) of initial thickness 1 mm in heated air (T = 1300° K) in the absence 
of straining. For this relatively thick fuel strip, ignition occurs at two distinct locations 
in an independent manner. Burning proceeds by each flame in an independent manner as 
well, evidenced by movement of each flame toward the adjacent oxidizer, in accordance with 
the stoichiometric requirements of the reaction. This ignition mode is consistent with the 
asymptotic analysis of ignition in [21] whereby relatively thick fuel strips at zero or low 
rates of strain ignite independently. Figure 2 demostrates that the interfaces burn in an 
independent fashion as well prior to the beginning of fuel burnout. 

Figure 3. in constrast, shows the evolution of temperature field associated with a fuel 
strip of the same initial thickness as in Figure 2 but exposed to a constant strain rate of 500 
sec -1 . While only a slight delay in ignition is observed in comparison to that in Figure 2, 
there is a very clear interdependence of the flames that do ignite; the flames then proceed 
to move toward each other, rapidly depleting the intervening fuel and extinguishing at a 
short time later. This behavior is a quantitative confirmation of the enhancement of thermal 
feedback suggested by the single step asymptotic solution of this ignition problem[21] and 
the corresponding dependent ignition mode for the fuel-air interfaces. The rapid burnout 
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process is further demonstrated by plots of propane, OH. and H radical mass fractions as a 
function of space and time shown in Figure 4. The figure also demostrates that the runaway 
in the H radical can precede the runaway in temperature during ignition of the interfaces, 
similar to the observations of Kreutz, et al. [9] for single strained fuel-air interfaces. 

Figures 5 and 6 show plots of temperature, propane, OH, and H radical mass fractions as 
a function of space and time for the propane fuel strip of the same initial thickness as above 
but with successively higher constant rates of strain. In Figure 5 there is a greater ignition 
delay than in Figure 4, but with ignition at a single location, indicating that ignition and 
sustained burning occur in the premixed or partially premixed mode. Here the rate of strain 
is so high that the fuel-air interfaces are brought together very rapidly, and substantial 
interchange of mass occurs before the ignition process. This mode of ignition again is a 
quantitative confirmation of the single ignition mode observed in the single-step asymptotic 
solution of this problem[21]. As is expected, this flame consumes the existing fuel-air mixture 
and burns out very rapidly. Figure 6, at a constant strain rate of 1500 sec -1 , demonstrates 
that at extremely high rates of strain, ignition can be prevented in the strained fuel strip; 
no significant formation of radicals, rise in temperature, or fuel consumption occurs. This 
“non-ignition” mode occurs due to the very rapid mixing of fuel and air, so high that fuel is 
diluted beyond the level required for ignition of the mixture. This observation is also seen 
in the asymptotic solution of the problem[21], and differs from modes of ignition prevention 
seen for single fuel-air interfaces in which heat loss due to high strain prevents thermal 
runaway. 

Finally, a sample result for solution of the viscous flow problem, equation (1), is shown 
in Figure 7a, with comparison made to the equivalent inviscid problem, equation (5), in 
Figure 7b. While the actual process of ignition and ignition times are the same for the two 
computations, the process of sustained burning of the fuel strip and eventual burnout is 
substantially longer for the viscous simulation than for the inviscid simulation. Clearly, the 
effect of viscosity is to reduce the effective velocity and strain field experienced by the species 
interfaces, transporting them more slowly toward each other and allowing for a reduced rate 
of consumption of intervening fuel. 


5 Summary and Conclusions 

The present study demonstrates that a full numerical simulation of ignition of a strained 
fuel strip with complete combustion chemistry exhibits the same phenomenological modes of 
thermal/radical runaway that are seen in a single-step asymptotic simulation of the problem. 
Further, burning and extinction processes may be explored via numerical simulation; these 
calculations suggest the sensitivity of the combustion characteristics to the degree of straining 
and the presence of viscosity in the flowfield. 
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Figure 1. Geometry of strained fuel strip problem, where ij a is the initial strip thickness. 
The strip is exposed to the temporally evolving strain rate e(t) in general. 
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Figure 3. Temporal and spatial evolution of temperature for the strained fuel strip with 
initial thickness 1.0 mm and constant rate of strain 500 sec -1 . The initial temperature of 
fuel is 400° K and of air is 1300° K. 



H Mass Fraction 


Fuel Consumption Rate 


0 oto r 


oooa \- 


0 006 f- 


0004 


000? k 


0000 



(EC 


*E-t 

46-4 

4W 

2E-* 

IE-7 

—1 OEO 


-0 50 0 00 0 50 

Distance (cm) 



Figure 4. Temporal and spatial evolution of temperature, fuel consumption rate, OH, and 
H mass fractions for the strained fuel strip with initial thickness 1.0 mm and constant rate 
of strain 500 sec -1 . The initial temperature of fuel is 400° K and of air is 1300° K. 
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Figure 5. Temporal and spatial evolution of temperature, fuel consumption rate, OH, and 
H mass fractions for the strained fuel strip with initial thickness 1.0 mm and constant rate 
of strain 1000 sec' 1 . The initial temperature of fuel is 400° K and of air is 1300° K. 
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Figure 6. Temporal and spatial evolution of temperature, fuel consumption rate, OH, and 
H mass fractions for the strained fuel strip with initial thickness 1.0 mm and constant rate 
of strain 1500 sec -1 . The initial temperature of fuel is 400° K and of air is 1300° K. 
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Figure 7: Temporal and spatial evolution of temperature for the strained fuel strip with 
initial thickness 1.0 mm and constant rate of strain 300 sec -1 . The initial temperature of 
fuel is 400° K and of air is 1300° K. 
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